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Abstract 

A study of multifragmentation of gold nuclei is reported at incident 
energies of 400, 600 and 1000 MeV/nucleon using microscopic theory. 
The present calculations are done within the framework of quantum 
molecular dynamics (QMD) model. The clusterization is performed 
with advanced sophisticated algorithm namely simulated annealing 
clusterization algorithm (SACA) along with conventional spatial cor- 
relation method. A quantitative comparison of mean multiplicity of 
intermediate mass fragments with experimental findings of ALADiN 
group gives excellent agreement showing the ability of SACA method 
to reproduce the fragment yields. It also emphasizes the importance 
of clustering criterion in describing the fragmentation process within 
semi-classical model. 
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1 Introduction 



A highly excited system formed in a nucleus-nucleus collision, as a rule, is 
expected to break into several pieces consisting of free nucleons, light charged 
particles (LCP's), intermediate mass fragments (IMF's) as well as heavier 
residue. This phenomenon of breaking of colliding nuclei into several pieces 
is known as multifragmentation [1, 2, 3, 4, 5, 6, 7]. Due to its complex 
dynamics, mechanism behind this picture of 'explosive' break up (into several 
entities) is not yet known completely. 

At low incident energies, excitation energy available to the system is very 
small. Therefore, larger impact of collisions is needed to break the system 
into pieces of different sizes. In other words, fruitful destruction is possible 
only for the central collisions. On the other hand, mutual correlations among 
nucleons are preserved in peripheral collisions, therefore, not much deviation 
from the initial picture will be seen. In contrast, excitation energy deposited 
in the system is very large at higher incident energies. Therefore, central 
collisions break the matter into very smaller pieces and rarely one sees inter- 
mediate mass fragments or heavy mass fragments in these events. Maximum 
number of IMF's can only be seen at semi-central impact parameters. Large 
number of experiments have witnessed this trend of fragmentation at various 
incident energies and impact parameters. This change in the behavior of 
fragment pattern is also termed as a rise and fall in the multifragmentation 
[6, 8, 9, 10]. 

As we go further towards higher incident energies, maximal of IMF multi- 
plicity starts shifting towards peripheral geometries. Such trends have been 
found and reported in several recent experiments of ALADiN - collabora- 
tion [2, 6]. In addition, manyfold aspects of spectator matter fragmentation 
have also been studied for the collision of ^^"^ Au -|-^^^ Au on ALADiN set up 
at incident energies varying between 150 and 600 MeV/nucleon. Recently, 
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INDRA experiments extended the energy domain covering the incident ener- 
gies between 40 and 150 MeV/nucleon [11]. The sole motivation for all these 
experiments was the fantastic physics that may emerge from the disintegra- 
tion of excited systems leading to the expansion of matter to low densities. 
This onset of multifragmentation and afterward transition to vaporization 
phase has also been linked to the concept of liquid-gas phase transition of 
nuclear matter [7, 12, 13]. Such critical behavior is, however, reported to be 
influenced by the finite size effects [3, 14] . 

All these experimental studies characterize the evolution of heavy-ion re- 
action from dominant multifragment-decay channel to complete disassembly 
into light charge particles (LCP's) and free nucleons sometimes also termed 
as 'vaporization'. The very recent study by Puri and Kumar [15] analyzed the 
^^Ca+^^Ca reaction for incident energies between 20 and 1000 MeV/nucleon 
and over entire impact parameter range. They predicted a clear rise and fall 
of multiplicity in incident energy and impact parameter planes. 

On the theoretical front, not much success has been reported to repro- 
duce the ALADiN experimental data [6, 16, 17]. Theoretical approaches 
which follow the evolution of target and projectile to complete disassembly 
of nuclear matter needs secondary algorithm to clusterize the phase space. 
Even afterburners have also been employed to extract fragments. The present 
study aims to check whether microscopic reaction models can explain the uni- 
versality reported by ALADiN group [8] in spectator fragmentation or not. 
Molecular dynamical models QMD [4] and QPD [18] were found to explain 
some of the features of this experimental data [6]. This questions the valid- 
ity of molecular dynamics models (MDM). The fallacy was largely attributed 
to the lack of advanced secondary clustering models [19, 20, 21]. The clus- 
tering criterion is one of the basic ingredients that may control the reaction 
mechanism in semi-classical models like quantum molecular dynamics model. 
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Recently, a novel clusterization algorithm based on the energy minimiza- 
tion criteria namely simulated annealing clusterization algorithm (SACA) 
was proposed [20]. As a first attempt, results with this algorithm were quite 
promising one [20]. In ref. [20], ^^"^Au +^^'' Au reaction was studied at inci- 
dent energy of 600 MeV/nucleon. Based on the ALADiN results, there one 
assumed that fragment pattern does not change above 400 MeV/nucleon. 
Therefore, it remains to be seen whether QMD model can reproduce this 
universality feature or not. We plan to address this situation in this letter. 
We apply this algorithm to ALADiN data at incident energies of 400, 600 
and 1000 MeV/nucleon in order to see whether our approach can explain 
the rise and fall phenomenon and universal behavior in spectator fragmen- 
tation at such higher incident energies. It is worth mentioning that SACA 
algorithm has been robust against experimental data at lower tail of incident 
energies. In our earlier studies [22] , SACA method was reported to reproduce 
the charge yields at incident energies between 25 and 200 AMeV. In this anal- 
ysis, 0-|-Ag/Br reactions were taken [22]. In another study, SACA method 
was tested against INDRA experimental data at 50 AMeV [23] . In this study, 
Xe-|-Sn reaction was subjected to multifragmentation and various variables 
such as charge, proton like and IMFs yields, angular distribution, average 
kinetic energies etc. were analyzed. SACA method explained all these ob- 
servables quite nicley, whereas conventional method failed badly [23]. Due to 
the fact that interaction energy among fragments is ignored, this approach 
of SACA can not be applied to incident energies below above mentioned one. 
To study fragmentation in Au+Au reaction, we followed nuclear collisions 
within QMD model [4]. The phase space thus generated is clusterized using 
advanced SACA method. 
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2 SACA formalism 



To generate the phase space of nucleons, we use quantum molecular dynamics 
(QMD) model. For the details of the QMD model, reader is referred to 
refs. [4, 20]. The next essential step is to clusterize the phase space stored at 
various time steps in each event. The extensively used approach assumes the 
correlating nucleons to belong to same fragment if their centers are closer than 
4 fm i.e. | — r/3 |< 4/m. It may often lead to wrong results if applied at 
higher densities and hence can't address the time scale of multifragmentation. 
This approach is labeled as minimum spanning tree (MST) algorithm. 

In our latest approach, fragments are constructed based on the energy 
correlations. The pre-clusters obtained with the MST method are subjected 
to a binding energy condition [20, 24]: 



with Ebind = -4.0 MeV if iV/ > 3 and Ebind — otherwise. In eq. (1), 
is the number of nucleons in a fragment and is the center-of-mass 
momentum of the fragment. The requirement of a minimum binding energy 
excludes the loosely bound fragments which will decay at later stage. 

To look for the most bound configuration (MBC), we start from a random 
configuration which is chosen by dividing whole system into few fragments. 
The energy of each cluster is calculated by summing over all the nucleons 
present in that cluster using eq. (1). Note that we neglect the interaction 
between the fragments. The total energy calculated in this way will differ 
from the total energy of the system [24]. 




'Bindi 



(1) 
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Let the total energy of a configuration k be Ek{= J^iNfCi)^ where Nf 
is the number of nucleons in a fragment and Q is the energy per nucleon 
of that fragment. Suppose a new configuration k' (which is obtained by 
(a) transferring a nucleon from randomly chosen fragment to another frag- 
ment or by (b) setting a nucleon free, or by (c) absorbing a free nucleon 
into a fragment) has a total energy E/^/. If the difference between the old 
and new configuration AE{— Eg — Ek) is negative, the new configuration 
is always accepted. If not, the new configuration k' may nevertheless be 
accepted with a probability of exp{—AE/v), where v is called the control 
parameter. This procedure is known as Metropolis algorithm. The control 
parameter is decreased in small steps. This algorithm will yield eventually 
the most bound configuration (MBC). Since this combination of a Metropolis 
algorithm with slowly decreasing control parameter v is known as simulated 
annealing, so our approach is dubbed as simulated annealing clusterization 
algorithm (SACA). For more details, we refer the reader to ref. [24]. 

3 Results 

For the present study, we use a soft equation of state (EoS) along with 
standard energy-dependent n-n cross section [25]. The soft EoS has been 
advocated by many studies [5, 6, 16, 17, 20, 26]. The phase space is generated 
and stored at many time steps and is then subjected to the above mentioned 
clusterization procedures. To address the time scale of multifragmentation of 
spectator matter, we employed SACA method as well as spatial correlation 
method {i.e. MST). 

The density of environment is often correlated with the prediction of 
breaking of nuclear matter into pieces. One can also look density distribu- 
tion in coordinate space to investigate the formation of fragments. We here 
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compute average density of system as : 

/ 1 At+Ap At+Ap 



1 




(27rL)3/2 



(2) 



with Yi and being the position coordinates of i and j nucleons. 
The Gaussian width L is fixed with standard value of 1.08 fm. Figure 1 (top 
panel) shows the time evolution of average nuclear density (p/ p^) for Au+Au 
system at incident energies of 400, 600 and 1000 MeV/nucleon and at an im- 
pact parameter of 6 fm. The average nuclear density reaches its maximal 
around 25 fm/c. This time domain also witnesses the maximum collision 
rate and nuclear interactions which are going on between among target and 
projectile nucleons. This maximal density shifts towards later times as we 
go down the incident energies. The fine point is that there is an insignificant 
change in the density profile while enhancing the incident energy by the fac- 
tor of 2.5 times i.e. going from 400 to 1000 MeV/nucleon. At the final stage 
of the reaction, we don't see any significant change with the incident energy. 
The bottom panel of fig. 1 shows the time evolution of the heaviest fragment 
^j^max\^ using MST and SAGA techniques. The MST method gives one big 
cluster at the time of maximum density, whereas one sees striking ability of 
SAGA method in identifying the heaviest fragment quite early when violent 
phase of the reaction still continues. This suggests that evolution of multi- 
fragmentation is an intricate process. In other words, fragmentation starts 
at quite early stage when nucleons are still interacting among themselves 
vigorously. The early recognition of heaviest fragment l^A^"^^) rules out its 
formation out of the neck region, i.e. geometrical overlap between projectile 
and target. This suggests the emission of l^A'"'"-'^) from the spectator region. 
Similar trends of transition from the participant to spectator fragmentation 
has also been observed and reported by ALADiN-coUaboration [8] . This find- 
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ing also confronts the common standpoint of thermal origin of fragments i. e. 
fragments are created after the thermalization sets in. Further after violent 
phase of reaction is over {i.e. after 60 fm/c), binding energy of all clusters 
in SACA method is greater than Esind^ the minimum binding energy needed 
to bind the group of nucleons into cluster. Fragments after time 60 fm/c 
leave the reaction zone without nucleon-nucleon correlations being destroyed 
further. Hence fragment configuration obtained at the earlier time can be 
compared with experimental data. Strikingly, earlier detection of fragments 
(not shown here) at all incident energies upto 1000 MeV/nucleon gives us 
possibility to look into the n-n interactions when nuclear matter is still hot 
and dense. Further, one is also free from the problem of stability of frag- 
ments. The failure of MST method to detect the fragments also questions 
its validity at incident energies as high as 1000 MeV/nucleon. Simple cor- 
relations method fails to detect the fragments even at these high excitation 
energies. The further rise in (A"*"^) after 60 fm/c using SACA technique is 
due to the reabsorption of surrounding light fragments by the heavier frag- 
ments. We see that heavier (A™"^) survive at smaller incident energies than 
at higher incident energies. The capability of QMD model clubbed with 
SACA method is illustrated in fig. 2 where we display the mean multiplicity 
of intermediate mass fragments (Nimf) as a function of impact parameter 
of the reaction. Also shown are the results obtained with MST method. Our 
model calculations with SACA method are in close agreement with ALADiN 
data [8] for Am reaction at all incident energies 400, 600 and 1000 

MeV/nucleon. As seen in the fig. 2, we also achieved a reasonable reproduc- 
tion of the shape of impact parameter dependence of (Njmf)- Due to shallow 
minima sometimes, we also see second minima before 60 fm/c in peripheral 
collisions. We show also the calculations at these minima marked as (*). We 
see that fragment structure at these minima is further closer to the data. 
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Further the peak value of (Nimf) and the corresponding impact parameter 
b is also well estimated with QMD + SACA method. The prominent feature 
of the spectator decay is the invariant nature of the IMF distribution with re- 
spect to the bombarding energy. The SACA method successfully reproduced 
the universal nature of spectator fragmentation at all the three bombarding 
energies. It is worth interesting to note that these universal features observed 
in multifragmentation of gold nuclei persist upto much higher bombarding 
energies than explored in this work [27] . In contrary, the normal spatial cor- 
relation method fails badly to explain the production of intermediate mass 
fragments at all incident energies. This questions the vahdity of MST method 
in explaining the fragmentation pattern in heavy-ion collisions. 

4 Summary 

We have studied multifragment-emission in ^^'^Au +^^'^ Au reaction at inci- 
dent energies of 400, 600 and 1000 MeV/nucleon where ALADiN experiments 

showed the universality in the production of intermediate mass fragments. 
For this study, we employed QMD model clubbed with energy minimiza- 
tion algorithm (SACA) along with conventional spatial correlation method. 
Our findings reveal that SACA is able to reproduce the universal nature of 
multifragmentation of excited spectator over entire impact parameter-energy 
plane whereas spatial correlation method failed to reproduce the IMF mul- 
tiplicity. This is for the first time that QMD + SACA approach is able to 
reproduce the entire energy domain. It also shows that mass and multiplic- 
ity of spectator fragments remain invariant to range of bombarding energies. 
This also resolved the earlier discrepancy where QMD model underestimated 
the fragment yield [6, 16] at large impact parameters even after 200 fm/c. In 
our case, SACA method is successful in breaking the spectator matter into 
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intermediate mass fragments. Our results show that the QMD model con- 
tains necessary ingredients to describe the physics of spectator decay. The 
clustering algorithm one uses, however, holds the key tenet to explain the 
reaction mechanism. 
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Figure Captions 

FIG. 1. Top panel: Time evolution of the average nucleon density (p/po) 
reached in ^^"^Au +^^^ Au collision. Bottom panel: The heaviest fragment 
(A""^) obtained with SAGA and MST analysis as a function of time in 
^^'^Au +^^'^ Au collision. 

FIG. 2. The mean multiplicity of intermediate mass fragments {Njmf) as 
a function of impact parameter b for the reaction of ^^'^ Au+^^'^ Au. The model 

calculations with SAGA (solid squares) and MST (open triangles) methods 
are compared with experimental data (open circles) reported by ALADiN 
group [8]. 
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